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SECTION I 


INTRODUCTION 

In a 1984 study the Energy Conversion Branch of the U.S. Air Force 
Wright Aeronautical Laboratories Aero Propulsion Laboratory investigated the 
use of advanced power generating devices for future mobile electric power 
(MEP) applications. One conclusion of that study was that kinematic 
Stirling engines have a high potential for meeting the requirements for mid- 
sized (30-100 kW) flightline and operational system electronics support 
applications. 

There are several kinematic Stirling engines currently under develop- 
ment by various organizations. An example is the automotive Stirling engine 
now being developed under the Department of Energy (DOE) Automotive Heat 
Engine Program with program management by NASA Lewis Research Center. One 
current effort also sponsored by this DOE program is Contract DEN3-351 with 
Stirling Thermal Motors, Inc. (STM) of Ann Arbor, Michigan, to experi- 
mentally evaluate advanced Stirling concepts [1]. STM, using primarily 
private funding, has designed, and is fabricating, a Base Technology 
Stirling Engine, designated the STM4-120 (4 cylinders, 120 cc swept volume 
per cylinder), incorporating these advanced concepts. The STM4-120 design 
is suitable for a variety of applications and since its heat input is from a 
liquid metal heat pipe it can be connected easily to almost any heat source. 

The STM4-120 design features several advanced concepts which should 
reduce the Stirling engine size, weight, complexity and expected manufac- 
turing costs. Indirect heating technology, incorporating heat pipes, is an 
integral part of the STM4-120, making it possible to divide the engine into 
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an energy conversion unit (ECU) and a distinctly separate externa! heating 
system (EHS), thus permitting simplification of the heat exchanger design. 

The heater of the ECU, actually the uniform temperature condenser of the EHS 
heat pipe, is designed to take advantage of the high film coefficient of the 
condensing metal vapor. This permits design optimization based on Stirling 
cycle thermodynamic requirements, and not on flue gas heat transfer, 
resulting in improved engine performance. 

Other advanced features, such as the variable angle swashplate power 
control and compliantly mounted reciprocating seals, are presented in the 
literature [2]. 

An analysis by Argonne National Laboratory looked specifically at the 
STM4-120's Air Force MEP applicability and concluded that such use is 
attractive [3]. 

These two studies and the STM4-120 design, led the Air Force Wright 
Aeronautical Laboratories to sponsor an additional task to the current DOE 
sponsored activity at STM. The new Air Force sponsored task was to design 
the external heating system for the STM4-120, which would meet the appro- 
priate military specifications, MIL-5TD-633 and MIL-G-52884 for a 30 kW 
mobile electric power generator set. 

The STM4-120 was designed to produce about 40 kW (53 HP) at 3000 RPM. 
The Air Force has an MEP requirement for a range of 20 to 60 kW. Although 
the STM4-120 will produce approximately 25 kW (33 HP) at 1800 RPM, the 
heating system will be designed to the full rated power of 40 kW. The 
eventual reoptimization of the STM4-120 to produce the military required 
power at 1800 RPM is a straightforward procedure and will have relatively 
little impact on the external heating system design. A complete description 
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SECTION II 

EXTERNAL HEATING SYSTEM CONCEPTUAL DESIGN 

1. EXTERNAL HEATING SYSTEM CONFIGURATION 

The conceptual design of the external heating system (EHS) was based on 
the heat input requirements of the STM4-I20 Base Technology Stirling Engine. 
The EHS was designed to operate on liquid fuel, using a flat plate, counter 
flow recuperator and a corrugated, flat-surface heat exchanger. In both 
components the flow of air and flue gas is strictly laminar (Poiseuille flow 
between flat plates) and hence the Nusselt number is independent of the 
Reynolds number. The Nusselt number based on gap hydraulic diameter, is 
taken to be always 7.6 [4]. 

The combustion chamber design was chosen from an existing Philips 
Stirling engine, the ADVENCO 4-88. The burner, atomizing air nozzle and 
mixing chamber were readily adaptable to the EHS design concept regarding 
fuel, size and performance. 

In operation, preheated air is mixed with atomized fuel in the combus- 
tion chamber and burns continuously. The atomized fuel is introduced into 
the combustion chamber through a nozzle. The geometry of the combustion 
chamber was designed for proper mixing and distribution of flue gas flow 
into the heat exchanger. Primary air is fed into the combustion chamber 
through a swirl chamber and secondary air is fed downstream through radial 
holes. Film cooling further downstream prevents overheating of the 
combustion chamber walls. 

Hot flue gas flows from the combustion chamber through the multiple 
parallel gaps of the heat exchanger, transferring its heat through the wall 
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into the finned enclosure. The enclosure is the evaporator section of a 
heat pipe where liquid sodium in the wick lining the inside wall evaporates. 

The sodium vapor flows to the engine via connecting tubes. 

Since the evaporator enclosure was designed to operate at sub-atmospheric 
pressure, coarse wire mesh was incorporated to provide structural support to 
the fin walls, while providing a vapor path for the sodium working fluid. 

The residual heat of the flue gas leaving the heat exchanger is used to 
preheat the incoming combustion air in the recuperator. 

A schematic diagram of the EHS is shown in Figure 1. Due to the use of 
heat pipes for heat transport the EHS can be located a variety of locations 
relative to the STM4-120 ECU. 

An analysis and simulation code for the EHS was developed to determine 
the influence of various geometric parameters on the system performance. 

The range of parameter variation was limited by the following 
considerations: 

• Packaging requirements limited the overall system dimensions to 
about 330 mm x 250 mm x 100 mm as shown in Figure 2. 

• Heat exchanger (evaporator) fins and preheater plates were 
designed with the same depth for ease of fabrication and packaging 
convenience. 

• In laminar flow heat exchangers the Nusselt number can be assumed 
to be independent of the Reynolds number, therefore it was 
advantageous to divide the flow path into as many parallel paths 

as possible, to derive the maximum reduction in flow friction 
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power loss, as shown in Table 1. The narrowest gap deemed 


practical was 1.0 mm. 


2. SYSTEM ANALYSIS AND OPTIMIZATION 

The geometrical parameters of the recuperator (or preheater) and heat 
exchanger alone determine the performance of the EHS provided that proper 
mixing and residence time can be achieved in the combustion chamber. Since 
an existing design of proven satisfactory performance was specified for the 



burner, the combustion chamber design was decoupled from the system analysis 
and optimization. 

a. Description of the System Analysis Code 

The purpose of the analysis and simulation code was to determine 
the EHS efficiency and the effect on flow losses and cold start fuel 
consumption based on the geometric parameters of the recuperator and heat 
exchanger. 

Aerodynamic and heat transfer models were employed to define the design 
parameters based upon the heat transfer between the exhaust and intake air 
in the recuperator, the combustion characteristics of the fuel and burner, 
and the heat transfer from the flue gas to the heat pipe evaporator [5] [6]. 



In heat transfer calculations, the temperature dependent properties of 
the air and flue gas were taken into consideration [7] and therefore, a 
finite differences numerical scheme was used to solve the energy equations. 
Furthermore, in calculation of the flame temperature and heat exchanger 
performance, the dissociation of the flue gas was taken into consideration. 











This effect was neglected in the analysis of the recuperator since the 
temperatures are sufficiently low. 

Since the total laminar pressure drop in the heat exchangers was 
expected to be less than 50 cm of water it was deemed unnecessary to take 
into account the variable fluid properties in the flow direction. The 
pressure drop was determined for each component separately, based on 
Poiseuille flow at the average temperature prevalent in that component. 
These average temperatures were calculated by numerical integration. 

The code was designed to accommodate any fuel (atomic hydrogen to 
carbon ratios between 0.25 and 4.0) and equivalence ratios ranging from 1.0 
to 1.4. 

The governing equations for the preheater and heat exchanger/evaporator 
heat transfer analyses are included in Appendix B. A sample computer 
run is included in Appendix C. 


b. Results of the System Analysis and Optimization 


In order to investigate the influence of the ratio of the evapora- 
tor area to the recuperator area (with a constant sum thereof) a series of 
simulations was made with the following geometric parameters unchanged; 
fin/plate depth, number of gaps and gap width in both the evaporator and 
recuperator, and the sum of the lengths of the evaporator and recuperator. 

The results are summarized in Figures 3 and 4 showing, respectively, 
the effect of the area ratio on the efficiency and the cold start fuel 
consumption. The area ratio is defined as the ratio of recuperator area to 
heat exchanger/evaporator area. The cold start fuel consumption decreases 
with an increase of the area ratio, and efficiency exhibits a maximum at an 
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FiRure 3 


Effect of Area Ratio on EHS Efficiency 
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area ratio value of approximately 0.25. A design point of 0.32 was chosen 
to account for possible deviation of the analytical model from reality and 
to allow sufficient interface space. That area ratio was held constant for 
the rest of the investigation. 

The effect of total heat transfer area at a constant ratio of heat 
transfer area to cross sectional area was investigated by a series of 
simulations with increasing depth of the evaporator fins and recuperator 
plates maintaining all other geometric parameters constant. The results are 
summarized in Figures 5, 6, and 7, showing respectively, the effect of the 
total heat transfer area on the efficiency, flow losses and cold start fuel 
consumption. A design point at the total area of 4,235 m corresponding to 
a fin/plate depth of 70 mm was selected as a favorable compromise between 
efficiency and flow losses on one hand and cold start fuel consumption 
(size, weight, and cost) on the other hand. 

It should be noted that the efficiency and cold start fuel consumption 
in Figures 3 and 4 differ from the values in Figures 5 and 7. This is due 
to the iterative procedure used to determine first the design point with 
regard to area ratio, then the design point with regard to total heat 
transfer area. The final value of total area is larger than the value used 
to determine area ratio. 

Finally, the sensitivity of the performance to the gap width was 
investigated by a series of simulations with varying gap width. The number 
of gaps for each simulation was determined based on the gap width so that 
the total width of the heat exchangers was maintained constant at 250 mm. 
The results of this investigation, summarized in Figures 8 and 9, shows as 
expected, that increasing the gap width has an adverse effect on the 
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efficiency which is much more severe than the beneficial effect on flow 
losses. Consequently a gap width of 1.0 mm was retained as the design 
point. 

The entire investigation was performed at full load corresponding to 
DF2 fuel consumption of 2.67 grams/sec with an equivalence ratio of 1.25 
(air /fuel ratio of 18:1), 4% atomizing air and no EGR. The material wall 
thickness was kept at 0.2 mm. The goals of this investigation were to 
establish preliminary design parameters, hence the effects of fouling and 
corrosion were not taken into consideration. 

The performance at the above conditions with the selected geometric 
parameters is summarized in Table 2 and Figures 10a and 10b. 
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Table 2 


STM4-120 EHS Performance Data 


Preheater 


Evaporator 


£ 

; 

£ 


Geometric Data 


Air Side 


Flue Side 


Length, mm 
Depth, mm 
Gap width, mm 
Number of gaps 
Wall thickness, mm 


Temperature data 

Average (bulk) temperatures, K 


220 

220 

110 

70 

70 

70 

1.0 

1.0 

1.0 

104 

104 

67 

0.2 

0.2 

0.2 

694 

827 

1370 


Fuel Data (DF2) 

Atomic H/C ratio 
Specific lower heating value 


1.804 

10.179 cal./gr. 


Operating Data 

Equivalence ratio 
Atomizing air 
Fuel consumption 
Ambient temperature 
Heat pipe temperature 


1.25 

4 * 

2,67 gr./sec. 
301 K 
1085 K 


Performance 

Preheated air temperature 

Flame temperature 

Heat exchanger outlet temperature 

Exhaust outlet temperature 

Heat output 

EHS Efficiency 

Cold start penalty 

Pressure drop 

Flow losses 


1034 K 
2397 K 
1097 K 
498 K 
104 kW 
91.6% 
128.3 gr. 
4.64 k Pa 
819 Watt 
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HEAT PIPE PRELIMINARY DESIGN 

The conceptual Stirling powered electrical generator set includes the 
STM4-120 Energy Conversion Unit (ECU) coupled to the EHS through four heat 
pipe tubes with mechanical couplings that simultaneously supply the ECU with 
heat via sodium vapor and return the sodium condensate to the EHS. 

In this section, a brief review of current heat pipe theory is 
presented, followed by a description of the EHS heat pipe configuration and 
the application of the theory to the particular design for the EHS. 

1. Heat Pipe Theory 

The theory of heat pipes is well documented in the literature, hence 
only a brief review of current theory is presented to provide a basis for 
the hydrodynamic and thermal design of the EHS heat pipe. Parts of this 
section are excerpted from "Theory and Design of Variable Conductance Heat 
Pipes," by D.B. Marcus, NASA Report Number CR-2018. A listing of other 
pertinent references is included in that report. 

a. Hydrodynamics 

Because a heat pipe involves the circulation of a working fluid, 
certain pressure drops arise. In general, there will be viscous losses due 
to liquid flow in the wick or capillary structure, and viscous and inertial 
losses due to vapor flow in the core. In addition, there may be body forces 
that either aid or hinder circulation (e.g. acceleration fields due to 

22 










etc.). For steady state operation of a heat pipe, a pressure head equal to the 
sum of these iosses must be supplied by the capillary in the wick. This yields 
the following steady state pressure balance, which must be satisfied between 
all points along the heat pipe. 


AP c = AP^ + AP y ± AP b (1) 


net I 

liquid 

vapor 


body force 

capillary = 

pressure + 

pressure 

± 

head 

. head J 

- drop 

. drop . 


. (if any) _ 


(1) Capillary Head 

Of the terms in Eq. 1, the liquid and vapor pressure drops are 
functions of the circulation rate and increase with the heat transfer load on 
the device, while the body force term is usually independent of the load. Thus, 
to satisfy Eq. 1, the capillary head must also increase with load in such a way 
as to match the losses incurred. 

The capillary head in a saturated heat pipe wick arises as a dynamic 
phenomenon. It is due to the existence of a pressure difference across a curved 
liquid-vapor interface which is given by: 


APi 


= 0 (- 


R. 


R. 


( 2 ) 


where APi = Interfacial pressure difference 
o = Surface Tension 

Rj, R 2 = Two orthogonal radii of curvature of the interface 



This pressure difference, which for concave menisci results in a depres- 
sion of the liquid pressure with respect to the vapor, exists all along the 
heat pipe wick. In order to obtain a net capillary head, it is necessary for it 
to be greater at the evaporator than at the condenser. In a heat pipe under 
load, this is exactly what occurs due to changes in the interface curvatures. 
Vaporization of the liquid in the evaporator causes the menisci to recede into 
the wick resulting in a decrease in the radii of curvature, while condensation 
in the condenser has the opposite effect. Therefore, capillary pressure is not 
constant. Capillary pumping is a passive phenomenon which automatically adjusts 
to meet the flow requirements, within limits. 

The extremes of the interfacial pressure difference, which can be gener- 
ated along the wick depend on the wick geometry and the wetting angle between 
the liquid and wick material. For the EHS application, the wick was specified 
to be stainless steel wire mesh, based on previous experience. For this type of 
wick, the maximum capillary pressure becomes: 


> = 
c max r c 


(3) 


w + d 


Where r c = effective capillary radius, = -2— ^ 
W = wire spacing 
d = wire diameter 


(2) Liquid Pressure Drop 


Liquid flow in heat pipes is generally characterized by very 
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low velocities ar.-d Reynolds numbers. Consequently, inertial effects can be 
neglected for steady state operation and the flow losses attributed only to 
viscous shear. 


The pressure drop in wire mesh wicks can be expressed as a form of Darcy's 
Law for flow in porous media, which is usually expressed in terms of a permea- 
bility - K, a measure of the wick flow resistance. 



In Eq. 5, d is the wire diameter and e can be calculated by the 
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equation: 


e = 1 - w/4 SNd 


( 6 ) 


where N = mesh number 

S = crimping factor «1.05 


(3) Vapor Pressure Drop 


The vapor pressure drop in heat pipes is often considerably 
more difficult to calculate than that in the liquid, for in addition to 
viscous shear, the analysis must account for momentum effects and perhaps 
turbulent flow and compressibility. Complicating this analysis is the fact 
that mass addition in the evaporator and mass removal in the condenser can 
significantly alter the velocity profiles and hence, the local pressure 
gradient. 

From conservation of momentum, the vapor pressure drop can be expressed 
as follows: 


«Py -»y Re y> My n»y 


2m. 


dm. 


dx 


2 \ V P v 


-e 


A 2 p 
v y v 


dx 


(7) 


Here, the first group represents the pressure drop due to viscous shear 
and the second group represents the pressure drop due to inertial effects, 
where f is the frictional drag coefficient and 0 is defined by the 
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Equation 8 accounts {or the effects of changes in vapor velocity across 
the section of consideration. 

(4) Body Force Head 

The last term in the pressure balance equation is the 
pressure head due to body forces acting on the working fluid. In the EHS 
heat pipe, the only body force that arises is due to the acceleration of 
gravity. Since the density of the vapor is very much lower than that of the 
liquid, the body forces on the vapor can be neglected. The pressure head 
due to the gravitational field is expressed by the equation: 



AP =J dP g =±J Po g cosOdz (9) 
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where AP g = cemponen, of the hydrostatic head in the 
direction of liquid flow 
g = acceleration of gravity 

0 = angle of the liquid flow axis with respect to 
the direction of gravitational acceleration 
dz = elemental length 
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(5) Capillary Pumping Limit 


The fact that there exists a maximum capillary head for any 
wick-fluid combination results in a hydrodynamic limit on heat pipe 
capacity. As mentioned previously, the capillary head must increase with 
the liquid and vapor pressure drops as the heat load increases. Since there 
exists a limit on the capillary head there must also exist a corresponding 
limit on the heat load if the pressure balance criterion is to be satisfied. 
This defines the capillary pumping limit, expressed by the following 
inequality: 


AP - AP S AP 0 + AP (10) 

c max g a v 

(6) Entrainment Limit 

A phenomenon which can effect the capillary pumping limit, 
but was not included in Eq. 10, is liquid entrainment in the vapor. In a 
heat pipe, the vapor and liquid generally flow in different, often opposite, 
directions. Since they are in contact at the wick surface, this sets up a 
mutual drag at the interface. If the relative velocity between the liquid 
and the vapor becomes too great, the interface becomes unstable and droplets 
of liquid will be torn from the wick and entrained in the vapor. Since this 
liquid never reaches the evaporator, it cannot contribute to the heat 
transferred by the heat pipe. However, it does contribute to the liquid 
flow loss. Thus, the maximum axial heat transfer in the heat pipe is no 
longer equal to the maximum fluid circulation rate times the latent heat of 
vaporization, but some lower value which defines the entrainment limit. 
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The conditions leading to entrainment are expressed in terms of the 
ratio of vapor inertial forces to liquid surface tension forces, called the 


& 


Weber number: 


P v V 2 z 


where p = vapor density 


= average vapor velocity 


liquid surface tension 


= characteristic dimension associated with liquid 


surface 


Limited experimental data with screen wicks indicate that a Weber 


number of unity represents the entrainment condition when the charac- 


teristic dimension, z, is set approximately to the screen wire 


diameter [8]. 


When the Weber number is set equal to unity, the limiting axial 


heat flux corresponding to incipient entrainment is given by: 


Py 

*■ z 


where Q = Maximum heat flux for incipient entrainment 

Cir»av 


A v = Vapor core flow area 


X = Latent heat of vaporization 





(7) Sonic Limitation 


It can be shown analytically [9] that there is a corres- 
pondence between constant area flow in a heat pipe with mass addition 
(evaporator) and removal (condenser) and constant mass flow in a converging- 
diverging nozzle. The end of the evaporator corresponds to the throat of 
the nozzle. Consequently, just as there is a sonic limitation on the flow 
velocity through a nozzle throat, there is a similar limit on the flow 
velocity at the heat pipe evaporator exit. For a given exit temperature and 
working fluid, this choked flow condition is a fundamental limitation on the 
axial heat transfer flux capacity of the heat pipe. 

The sonic limit is calculated by setting the vapor flow velocity equal 
to the sonic velocity in the continuity equation and multiplying by the 
latent heat of vaporization as follows: 


max 


= A p V 
v s 


(13) 


where 


Os = axial heat flux at Mach 1 conditions 

max 

V s - sonic velocity of the vapor 


b. Heat Transfer 

The fluid circulation phenomenon discussed in the previous 
sections arise as a result of heat transfer into the heat pipe at the 
evaporator and out of it at the condenser. In most heat pipes this heat 
must be transferred through the walls and saturated wick. Generally, these 



processes are the major source of temperature drop in the heat pipe. 
Temperature drops also arise due to vapor flow losses along the heat pipe 
and to non-equilibrium at the liquid-vapor interfaces. 

(1) Boiling Limitation 

The fact that heat is transferred into the heat pipe through 
the wick gives rise to another limit on heat pipe capacity. Frequently, 
liquid is vaporized only at the wick surface as a result of heat conducted 
through the wick. However, the vapor at the wick surface is thermo- 
dynamically saturated, hence the fluid within the wick is superheated by 
virtue of the curvature of the menisci and the temperature gradient in the 
saturated wick. The greatest superheat occurs at the interface of the wick 
and wall. If this superheat becomes too large (it increases with the heat 
transfer rate), the fluid will begin to boil within the wick. This results 
in evaporator dryout due to a reduction in liquid flow area in the wick 
caused by nucleation. 

Further discussion of boiling limitation is referred to the literature, 
since the heat flux was specified to be a value experimentally proven to be 
within the boiling limitation of the EHS evaporator configuration. 

(2) Heat Pipe Temperature Characteristics 

Generally, heat pipe temperature characteristics are 
evaluated at full load, since the associated gradients are at a maximum 
under these conditions. The temperature drops associated with heat pipes 
occur in the evaporator wall, the saturated wick, and the condenser wall due 
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to conduction. There is an additional, generally much smaller, temperature 
drop associated with the vapor pressure drop between the evaporator and the 
condenser due to vapor flow resistance. 

(a) Temperature Drop Across the Evaporator Wall 

The temperature of the outer wall of the evaporator is 
essentially constant at 1085 K (812°C). The temperature drop due to conduction 
across the wall can be expressed by the following equation: 

_Q_ _ u dj 

"A " k dx (14) 


(b) Temperature Drop Across the Saturated Wick 

The temperature drop across the saturated wick can be 
evaluated using Eq. 14. However, the thermal conductivity must be replaced 
by an effective conductivity of the wick/liquid combination. The effective 
conductivity, k for saturated screen wicks is: 


M (k t * V - <‘ -eXKi-kJ] 
'« = l< k t * k w > ♦ (1 - eXk t - k w >] 

where k ^ = Liquid Thermal Conductivity 
k w = Wick Thermal Conductivity 
e = Wick Porosity 


m l o v vv ‘.*w^vvv.vvi w 


EHS. The heat pipe configuration is shown schematically in Figure 11. 

Liquid return is accomplished via gravity by locating the EHS heat exchanger 
below the ECU. 

Tests conducted at STM on the mechanical couplings, Aeroquip Conoseals, 
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Schematic Representation of STM4-120 Heat Pipe 
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available from the Marman Division of Aeroquip, proved the functionality and 
durability of these seals in sodium heat pipe applications. One heat pipe 
with this coupling has operated for more than 1775 thermal cycles of heating 
to 800°C and cooling to 150°C. Another heat pipe with its wick capillary 
connection across this coupling has operated against gravity for more than 
1000 hours. These tests are continuing [1]. An exploded view of a Conoseal 
joint is shown in Figure 12. 

The EHS heat pipe configuration was examined as both buffered (moving 
front) and non-buffered designs. The buffered configuration permits easy 
start-up and temperature control but suffers from the fact that the sharp 
temperature gradients in the front create extremely high thermal stresses, 
especially in the engine heater. As the front passed through this heat 
exchanger, some of the tubes would be at 800°C and others would still be at 
ambient, resulting in unacceptably high stresses, and increasing the 
probability of heater tube failure. 

In a non-buffered configuration, the entire heat pipe increases in 
temperature uniformly, eliminating the difficulties that result from thermal 
stresses. However, this configuration can experience start-up difficulties 
by condensing and freezing the sodium, drying out the evaporator before full 
fluid circulation can be established. 

The EHS heat pipe fluid circulation was designed to return the maximum 
amount of sodium to the evaporator upon shut-down by wicking only the 
evaporator. This technique is only applicable if liquid return from the 
condenser can be accomplished without capillary pumping. In the EHS heat 
pipe, liquid return is accomplished via gravity. 

Start-up is more readily accomplished since raising the temperature of 
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the working fluid occurs at the heat source, not at remote sections of the 
heat pipe which often require auxiliary heaters to melt the sodium to 
initiate fluid circulation. 

The liquid collects in the bottom of the evaporator fins, where it 
contacts the wick lining the inside walls of the fins. The wick supplies 
liquid sodium over the entire inside surface of the evaporator through 
capillary pumping. 

The evaporator configuration was based on the more difficult heat 
transfer from the flue gas to the evaporator wall, resulting in an average 
heat flux of 100 kW/m *. The maximum heat flux occurs at the bottom of the 
evaporator fins, where the flue gas enters the heat exchanger. The 
magnitude of the heat flux at this point is approximately 500 kW/m* . In 
experiments on sodium heat pipes at Philips and STM, this magnitude of heat 
flux was defined to be below the boiling limitation for sodium at 800°C, 
with commonly used wire mesh wicks. Once the configuration was established 
as compatible to the convective and radiative properties of the flue gas and 
the boiling properties of sodium, the design task was directed towards 
specification of a wick structure capable of maintaining an adequate supply 
of liquid sodium over the entire evaporator surface, i.e. satisfying the 
capillary pumping requirement, and providing structural support to the flat 
fins. 

A coarse wire mesh was placed in the vapor path to provide structural 
support for the fins at a slight penalty in temperature uniformity due to 
vapor flow resistance, while simultaneously providing positive contact 
between the capillary wick and the evaporator wall. A schematic cross- 
section of the evaporator is shown in Figure 13. 










3. Analysis of the Heat Exchanger/Evaporator 


An analysis was conducted on the evaporator to evaluate capillary 
pressure as a function of position in a typical fin using various wick and 
support structures, while maintaining the external geometry defined in Table 2. 


The analysis was developed from basic heat pipe theory and empirical 
relations for flow properties of wire mesh wicks as previously summarized 
[ 10 ], [ 113 , [ 12 ]. 

The most critical aspect of the evaporator design is maintaining suffix 
ciently high capillary pressure in the wick to distribute the liquid 
throughout the fins. In addition, the coarse supporting mesh must withstand 
the structural load due to the pressure difference across the fin walls 
without significantly restricting the vapor flow out of the fin. 

Sonic and entrainment limitations are a much lesser concern since the 


large number of fins, hence large vapor flow area results in vapor 
velocities well below a Mach Number of 0.2. This also simplifies the 
evaluation of the vapor pressure drop since the effects of the vapor 
compressiblity can be neglected. 

The heat flux from the flue gas to the evaporator varies according to 
the temperature distribution as shown in Figure 10. However, the evaporator 
analysis is greatly simplified by making the conservative assumption of 


constant heat flux. This assumption results in liquid flow only in the X 
(vertical) direction, and vapor flow only in the Y (horizontal) direction, 
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P c (x,y) 

P c (ref) 


= capillary pressure at a reference position (x ref , y ref ) 
= vapor pressure drop between (x,y) and (x fef , y fef ) 

= liquid pressure drop between (x,y) and (x^, y fef ) 

= hydrostatic head (due to gravitational body force) 



pressure must be equal to the sum of the pressure drops in the rest of the 
heat pipe. 

Vapor and liquid flow resistance are functions of the mesh geometry and 
fluid properties. Development of the liquid and vapor pressure drop 
equations 4 and 7, respectively, to the EHS finned evaporator application 
can be found in Appendix D. The results of the analysis are summarized in 
Figure 15. 
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satisfied the capillary pumping requirement, however, the chosen 
configuration was the best compromise among cost, manufacturability and 
performance. The results of the performance evaluation are summarized in 
Table 3. 
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Number of fins = 68 

Width/fin = 2.67 mm 

Wall thickness = 0.2 mm 

Wick structure = 4 layers of 165 mesh, 

.050 mm (.002”) wire diameter 

Support structure = 2 layers of 8 mesh, .48 mm (.017”) wire diameter 

Operating Data 

Heat pipe evaporator temperature 1085°K 

Heat input 104 kW 

Maximum wick capillary pressure 3200 Pa 


Performance 

Minimum capillary pressure 

650 Pa 

Maximum required capillary pressure 

2525 Pa 

Temperature drop due to: 

Conduction across evaporator wall/wick 

1.8 K 

Vapor flow restriction in fins 

0.4 K 

Vapor flow restriction in connecting tubes 

2.0 K 

Condenser (engine heater) 

4.8 K 

Total temperature drop 

9.0 K 

Temperature of inside of wall engine heater tubes 

1076°K 
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SECTION IV 

EXTERNAL HEATING SYSTEM CONSTRUCTION 

The construction of the EHS was designed to be adaptable to simple 
folded sheet metal manufacturing techniques. The folded constructions are 
amenable to robotic welding techniques. In this section, a brief 
description of the envisioned construction procedure is presented. 

1. Construction of the EHS Heat Exchanger/Evaporator 

The evaporator construction begins by bonding the layers of wire mesh 
forming the wick to the sheet that forms the outer wall. The bonded sheet 
is folded over mandrels of 1.0 mm thickness for the flue gas gap, and 2.0 mm 
for the sodium vapor gap, resulting in a structure like that shown in Figure 16. 
The folded sheet is then soaked in solvent to remove the adhesive and any 
oil or grease. 

The folded sheet is then mounted on a special mandrel to perform the 
squaring operation on the root of the fins as shown in Figure 17. In 
addition, the edges of the fins are crimped together. This operation 
permits ready welding of the fin edges without filler rod and simplifies the 
welding of the roots of the fins to the rest of the enclosure. 

The fins are filled with the coarse supporting mesh, then mounted in a 
welding fixture with heat sinks. At this stage, the fin seams are welded. 

Once the fin seams are welded, the assembly is placed in a different 
welding fixture and the top and bottom covers are welded in place. These 
covers fit accurately into the squared fin roots, again permitting welding 
without the use of a filler rod. See Figure 18. 
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The final step consists of filling the evaporator plenum, formed by 
the sides of the outer fins and the top and bottom covers, with coarse 
supporting mesh, then welding the back plate in place. The pack plate 
contains the four connecting tubes, with mechanical coupling flanges, that 
supply the STM4-I20 ECU with sodium vapor. 

2. Construction of the EHS Recuperator 

The techniques used to form the recuperator passages are identical to 
the initial folding, squaring and seam welding operations used to form the 
evaporator enclosure. 

3. EHS Preliminary Design Drawings 

The preliminary design drawings for the complete EHS are included in 
Appendix E. 

4. Heat Exchanger/Evaporator Scale Model Tests 

In order to confirm the functionality of the novel evaporator design, 
scale models were built using the construction techniques outlined above. A 
photograph of the model is shown in Figure 19. 
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SECTION V 


CONCLUSIONS 


The conceptual design of the STM4-120 External Heating System appears 
to be a relatively simple, efficient, and inexpensive system for supplying 
heat to the Stirling Energy Conversion Unit to provide shaft power for a 
mobile electric power generator set. In addition, the preliminary design 
and subsequent scale model testing of the EHS heat exchanger/evaporator 
substantiate the viability of the conceptual design. 

The system optimization resulted in an EHS package size of 330 mm x 
250 mm x 70 mm, operating at a full load output of 40 kW with an equivalence 
ratio of 1.25, 4% atomizing air, fuel consumption rate of 2.67 g/sec, and an 
EHS exhaust temperature of 225°C. 

The analysis of the heat exchanger/evaporator and scale model testing 
indicate the viability of the evaporator configuration for supplying heat at 
800°C to the STM4-120 Stirling engine. The analysis also provides the 
framework for design algorithms for future designs. 

The primary concern of the heat exchanger /evaporator preliminary design 
was to specify a wick structure capable of providing sufficient liquid 
distribution within the evaporator fins at an average heat flux of 100 kW/m 2 
and a maximum of 500 kW/m 2 . 

Results of the preliminary design yielded an evaporator with 68 fins, 
each with a wick structure composed of four layers of 165 mesh, .050 mm 
diameter wire. The fin support structure was composed of two layers of 8 
mesh, .43 mm diameter wire. The specified wick, capable of producing 3200 
Pa of capillary pressure, was required to produce only 2500 Pa to provide 
adequate liquid distribution. 
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While much of the military MEP specifications, MIL-STD-633 and 
MIL-G-52884 is not applicable to a radically new and different prime mover, 
such as the STM4-120 Stirling engine, none of the existing design 
specifications are in conflict with these requirements. 

The authors feel that the results of this and related investigations 
warrant further development of the STM4-120 for Air Force MEP applications. 
Potential difficulties not accounted for in the preliminary design are 
problems related to fouling and corrosion. 

Counter flow recuperative preheaters have been developed by various 
manufacturers for gas turbines. These have very similar requirements to 
those of the EHS recuperator. It is the intention of STM to make use of the 
know-how of the gas turbine manufacturers, particularly in the areas of 
fouling and oxidation corrosion of the recuperator and heat exchanger/ 
evaporator, and also in addressing weld or braze joint leakage. 

Furthermore, it is possible that proven, available recuperators, 
metallic or even ceramic, can be directly applicable to the EHS system. 

This will reduce the technical risk and is thus the preferred approach 
provided that the system performance does not have to be excessively 
compromised. 

Solutions to these problems require further research and experimental 
testing and were not within the scope of this effort. 

The results of this investigation add further support to the 
conclusions drawn in the Argonne study, "Base Technology Stirling Engine 
Military Applications Assessment," listed below: 

• Fuel flexibility (similar to that possible with gas turbine sets); 
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• Low specific weight (comparable with the best achievable with 
diesel sets); 

• Low noise (better than that achievable with silenced diesel sets); 

• Low IR - both because the majority of heat is emitted through the 
radiator at already low temperature and because the engine system 
is more efficient than current sets and consequently releases less 
heat; and 

• High efficiency (comparable to the best achievable by diesels in 
this power range. 

Stirling Thermal Motors, Inc. is currently engaged in the fabrication 
of five prototypes of the STM4-120 ECU to establish the predicted 
performance and reliability of this Stirling engine configuration. This 
program is expected to continue into 1986, during which prototypes of the 
EHS design are anticipated to be built and tested. 
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INTRODUCTION 

A new concept In Stirling engine technology Is embodied In the 
engine now being developed at Stirling Thermal Motors, Inc. In Ann Arbor, 
Michigan. This Is a versatile engine suitable for many different applica- 
tions and heat sources. 

The engine, rated at 40 kW at 2800 rpm. Is a four-cyl Inder , double- 
acting variable displacement Stirling engine with pressurized crankcase 
and rotating shaft seal. It Incorporates Indirect heating technology with 
a stacked heat exchanger configuration and a liquid metal heat pipe 
connected to a distinctly separate combustor or other heat source. It 
specifically emphasizes high efficiency over a wide range of operating 
conditions, long life, low manufacturing cost and low material cost. 

This article describes the new engine. Its design philosophy and 
approach. Its projected performance, and some of Its more attractive 
appl ications. 


BACKGROUND 


In 1072, Ford Motor Company obtained a worldwide exclusive license 
from N.V. Philips of the Netherlands for the Stirling engine, covering 
virtually all applications, including automotive. 

Under this license agreement, the Research Lab at N.V. Philips 
was to design and build four 175 hp engines, two of which would be 
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In 1976 the two Stirling powered Torinos and an older Philips 
Stirling bus equipped with a 4-cylinder rhombic drive Stirling engine [2] 
were successfully demonstrated for three days in Dearborn, Michigan. 

A few years later, in 1978, Ford terminated its Stirling engine 
activities to make manpower available for short-term technological 
problems. A year later, Philips also stopped work on the Stirling engine. 

Upon these events, Stirling Thermal Motors, Inc. (STM) was founded 
in the United States to continue the work done at Philips, so that the 
results of the years of research and development work done at Philips 
since its last license agreement in 1968 - which had resulted in a 
technical breakthrough - would not be lost. 

STM's main purpose, from its foundation, was to develop commercial 
Stirling engines. Philips Laboratories had only made laboratory models for 
use in research. The only engine made for a special purpose was the one 
for Ford. When this particular engine was made. Philips was confronted with 
the practical reality of designing and building a Stirling engine for the 
most complex application imaginable - an automotive engine. During this 
time it was discovered that some components of the engine might form 
obstacles in the way of commercialization because of their complexity 
and vulnerability. 

From 1974 on, real breakthroughs were made in avoiding these 
complexities. This made a more simple four -cylinder, double-acting Stirling 
engine possible. 

Unfortunately, by this time it was too late to incorporate these 
improvements into the Ford engine. The intent was to use these new 
developments in a second-generation Ford engine. This, though, was not 
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done before Ford dropped its Stirling engine program. The engines being 


built by earlier Philips licensees were based on designs older than the 


Ford engine. Their configurations had been frozen for several years. It 


was therefore impossible to utilize the new improvements. 


From the outset STM was convinced that the time was ripe for 


commercialization of the Stirling engine because all the ingredients 


for a simple^ inexpensive and reliable engine with a long service life 


were present. 


GENERAL APPROACH, BASIC APPLICATIONS 


STM's general approach is based on the conclusion that competition 


with existing internal combustion engines should be avoided, at least in 


the beginning. Rather, markets should be found where the IC engine cannot 


be used and where the use of the Stirling engine would be very economical. 


making use of the unique properties of the Stirling engine. Of the many 


possible applications, particular attention was given to the following 


three: 


t Prime mover for heat-driven heat pump 


• Solar energy conversion 


• Engine for generator sets. 


If the manufacturing cost of the engine could be sufficiently low. 


particularly in mass production, the market in these fields alone could be 
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vast, 

TECHNICAL APPROACH 

V 

The whole drive of STM is to commercialize the Stirling engine. This 

4* 

means that the engine must be simple, reliable, inexpensive, and that it 
should have a long service life. None of these requirements should have an 
adverse effect on the performance of the engine. 

More than three years of designing, discussions with suppliers and 
vendors, component testing and price calculations, led to the new version 
of the Stirling engine. Studies done for NASA have shown that this engine 
configuration is suitable for a whole range of power sizes up to 500 hp. 

Special emphasis was placed on the flexibility of the engine to 
adapt readily to a wide range of specific applications, duty cycles and 
heat sources. 

Consequently, indirect heating technology is an integral part of 
the development effort, making it possible to divide the engine into an 
energy conversion unit (ECU) and a distinctly separate external heating 
system (EHS). Different heat sources coupled to the same ECU will adapt 
the engine to different applications and enhance commercial introduction 
since most of the development complexity and cost is in the energy con- 
version unit. 

c 

A liquid metal heat pipe [3] is used to transport the heat from the 
heat source to the expansion heat exchanger of the thermal converter. 

So far, most of the development effort has concentrated on the ECU, 
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Figure 5 

at exchangers of the S T H4- 120. There is one per cylinder. These will later be electron-beam 
e heat exchanger stack. The tubes are curved, enclosed in a flexible cannister, and 
o end Dlates. 












In 1974, during the work on the automotive engine, a relatively simple, 
heavy duty construction was found to vary the power [5], [6], [7]. In this 
case the mean pressure of the engine stays the same, but the stroke of 
the pistons changes. This method of power control has the further advan- 
tage of high part load efficiency. Such a construction could be used only 
with a swashplate drive since the stroke of the pistons Is controlled by 
the angle of the swashplate. It was tested thoroughly In a test rig and 
applied In the Advenco engine, but the Advenco engine was never adequately 
tested. Philips eventually sold It to NASA, where further testing was done. 

Figure 7 Illustrates schematically the variable swashplate mechanism. 
The swashplate Is mounted on a part of the shaft which Is tilted an angle a 
from Jte main shaft axis. The swashplate Is mounted In such fashion that Its 
centerline makes an angle a with the tilted shaft axis and It can be rotated 
relative to and about the tilted shaft axis In order to change Its angle and, 
with this, the stroke of the pistons. 



figure 7 - Schematic of the variable swashplate mechanism, showing the principle of changing the angle of 
the swashplate, which it makes with the shaft from 0 to 2 a. The small plate A is fixed on the shaft S with 
an angle a , The engine swashplate (drawn here as a solid line) is perpendicular on the shaft S. This situation, 
a - a • 0, means that the strokes of the pistons are zero. When the engine swashplate B is turned 180° relative 
to the small plate A, the angle then becomes represented by the dotted lines. In this case, the strokes of the 
pistons are maximal. By turning B relative to A between 0* and 1B0* any angle of the swashplate between 0 and 2a 
can be obtained, so the strokes of the pistons can be changed from zero to maximum. C is a hydraulic vane 
mechanism, the housing of which can turn in one or the other direction depending on which chambers are pressur- 
ized and which are not. The oil is fed via two channels in the shaft S. The turning of the housing is trans- 
mitted to the engine via bevel gears t and F. 
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Practical aodel of a variable swashplate, shown in two positions. 


The torque applied by the actuator to the swashplate in order to 


maintain a certain angular position depends on whether it was rotated 


to such position in the positive or in the negative direction. Rotation 


in the negative direction requires less torque since the engine torque 


Itself acts in this case to Increase the swashplate angle. 


Figure 10 shows the actuator torque as a function of the swashplate 


angular position relative to the shaft. The curves labeled and 

M (i^) refer to that torque for the positive and negative direction of 


rotation respectively. The third curve, y (<M shows the corresponding 


swashplate angle. The curve M (^) reverses its sign within its range of 


definition. The point of sign reversal is an unstable control point to 


be avoided by narrowing the range of definition to exclude it. In the 
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Rotary actuator torque requirements. H + means relative rotation for larger stroke in the sane direction 
as the rotation of the engine. N. means relative rotation for larger stroke in the opposite direction 
from the rotation of the engine. 


Energy Conversion Unit, the angle a is 12.5°, yielding maximum theoretical 
swashplate angle of 2a ■ 25°, corresponding to 180° rotation. The maximum 
swashplate angle of interest is only 22 e , corresponding to a narrower range 
of definition (124°) in which the torque M_(4>) does not change sign. 

Loss of hydraulic power will result in the gas forces bringing the 
swashplate to a position perpendicular to the main shaft axis, reducing the 
piston stroke to zero - an automatic safety feature. 
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SEALS 


In a four-cylinder, double-acting Stirling engine there are two 
types of dynamic seals: 

e Dynamic seals as piston rings to divide the four cycles from 
each other, and 


e Dynamic reciprocating seals on the piston rods, in order to 
contain the high pressure working gas in the engine. These 
seals should also prevent oil penetration Into the cylinders 
from the lubricated drive. 


For the dry-running piston rings, a good solution is found in 
using a reinforced PTFE material. 

However, the different types of reciprocating seals for the piston 
rods are still not reliable. Philips developed the rolling diaphragm seal, 
but this was shown, in the Ford/Philips engine, to be vulnerable in 
non-laboratory environments. 

STM was able to avoid the gas containment function of these 
reciprocating seals entirely. 

The new power control, with its variable swashplate, made it 
possible to enclose the relatively small drive with a pressure hull and to 
use a commercially available rotating shaft seal. Preventing oil penetration 
into the cylinders is, in this case, much easier and has already been 
thoroughly tested in other engines. 



SPECIAL FEATURES 


Amenability to dynamic balancing and the ease of starting the 
engine are two additional features of the variable swashplate drive 
and power control elaborated upon in this section. 

DYNAMIC BALANCING is achieved by adjusting the swashplate moments 
of inertia to the reciprocating mass. This is done in a manner enforcing 
perfect dynamic balance at a certain swashplate angle within its range 
of variation. At different angles unbalanced moments will appear, but since 
perfect balance automatically occurs at zero angle, these will be very 
small. 

STARTING of the engine can be accomplished by heating up the 
expansion heat exchanger and the regenerator and then suddenly using the 
control pressure to increase the swashplate angle. This will cause 
the pistons to move in their normal way, causing the engine to 
Immediately develop sufficient power to perpetuate the motion. An 
accumulator fed by the hydraulic pressure pump will be used for 
that purpose. 

Obviously, such a simple procedure may be used only for such 
applications as solar conversion since no accessories are required for the 
combustion. In other cases only a very small starter motor is required 
to power the accessories needed for combustion, such as the air blower. 

When the engine has reached its correct temperature the accumulator 
pressure may be used to quickly increase the angle, having the engine 
self-start as mentioned above. 
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DESCRIPTION OF THE ENGINE 


A layout drawing of the Energy Conversion Unit is shown in Figure 11. 
This unit is distinguished by two ir.ujor features: 

• Variable swashplate drive and power control contained in a 
pressurized "crankcase" with a commercially available rotating 
shaft seal containing the working fluid and making it possible 
to avoid the reciprocating rod sealing problem; and 

• Indirect heating, featuring heat-exchanger-stack configuration 
and liquid metal heat pipe heat transport system. 

Following is a brief description of key components. 

The rotating shaft seal assembly (1) is bolted to the rear crankcase. 
The drive shaft (2) is splined to the main shaft to illow for slight 
deflection or misalignment without affecting the seal performance. The 
drive shaft is supported radially by one needle and one ball bearing, but 
since it passes through the pressurized crankcase, it must support 
the corresponding thrust load. The thrust bearing (3) is a heavy duty tapered 
roller bearing. 

The rotating shaft seal (4) is a commercially available mechanical 
face seal similar to those found in compressors and gearboxes. The seal is 
secured in the housing by an interference fit. A rotating mating ring (5) 
is mounted on the shaft. Both components have an extremely flat, smooth 
sealing face. When the seal is installed, these faces are pressed together 



1) Rotating shaft seal assembly 

2) Drive shaft 

3) Rotating shaft seal 

4) Hating ring 

5) Radial lip seal 

6) Thrust bearing 

7) Pressure hull 

8) Rear main bearing 

9) Rear crankcase 


Figure 11 


The STH4-120 


10) Rotary actuator 

11) Hain shaft 

12) Bevel gear 

13) Swashplate 

14) Front main bearing 

15) Oil pump module 

16) Front crankcase 

17) Cooler 

18) Regenerator 




19) Heater 

20) Hot connecting duct 

21) Piston assembly 

22) Cylinder/regenerator housing 

23) Hydraulic service assembly 

24) Power control valve 

25) Oil scraper/capseal assembly 

26) Crosshead 

27) Accumulator 









and held in contact by compressed springs in the seal, then surrounded by 
oil. Because of their closely controlled surface finish, virtually no 
oil can pass between the sealing faces. 

To minimize wear, the seal is hydraulically balanced to reduce contact 
forces due to crankcase pressure and subjected to a steady flow of oil 
to maintain cool operation. 

The radial lip seal (5) prevents the oil surrounding the rotating 
shaft seal from draining into the sump. 

The pressure hull (7) fits over the rotating shaft seal housing 
and is bolted to the front of the crankcase, covering the entire drive. 

The pressure hull serves no structural function other than containing the 
crankcase pressure. 

The front (16) and rear (9) crankcase castings form the major 
building blocks of the engine, housing the entire drive and power control. 

In addition, the front crankcase houses the cold thermodynamic system 
including the coolers, water passages and cold connecting duct. Both 
are cast from Reynolds 390 Aluminum for reduced weight and suitable running 
surfaces for the reciprocating components. 

The swashplate (13), combined with the main shaft (11) provide 
the mechanism for converting the reciprocating motion of the pistons into 
rotary motion. The swashplate is mounted on the shaft in a manner such 
that relative rotation to the shaft causes the swashplate angle to change. 

The front (14) and rear (8) main bearings provide support for the 

main shaft. The two bearings are spherical roller bearings capable of 
accomodating the shaft deflection caused by the swashplate loads. 



AV. 
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The rotary actuator (10) is a hydraulic vane motor comprising two 
diametrically opposed vanes attached to the shaft and two attached to the 
housing, forming two chambers. Rotation of the actuator housing relative 
to the shaft is accomplished by pressurizing one chamber and relieving the 
other. The relative rotation between the shaft and actuator housing is 
transmitted to the swashplate via the bevel gear (12) to change the 
swashplate angle and therefore the piston stroke. 

The oil pump module (15) houses the lubrication and hydraulic pumps 
which are directly driven by the matn shaft. Both pumps are high-efficiency 
geroters. Internal passages duct the oil to galleries to provide lubrication. 
High pressure oil is ducted to the externally mounted hydraulic control 
system and through annular passages in the shaft to the rotary actuator. 

The hydraulic service assembly (23) is mounted to the underside of 
the front crankcase for ready access. It consists of the oil sump, hydraulic 
accumulator and the power control valve. 

The accumulator (27) maintains a large volume of oil at hydraulic 
system pressure to prevent lags in response to power demands. The accumulator 
permits use of a smaller pump, thus reducing parasitic losses. 

The power control valve (24) is a solenoid activated proportional 
spool valve that controls the pressure difference across the rotary actuator 
and in turn, the angle of the swashplate and piston stroke. 

The four coolers (17) are conventional tubular heat exchangers mounted 
in the front crankcase. There is one cooler for each cycle. 

The high temperature components consist of identical assemblies for 
each cycle forming a heat exchanger stack configuration. These components 
are investment cast from the lower cost iron-based superalloy CRM-60. 


The engine heater (19) is a simple tube bundle heat exchanger consisting 


of relatively short and small diameter tubes enclosed in a flexible cannister 
and brazed to end plates. The heater is then electron beam welded to the 
hot connecting duct (20) and cyl i nder /regenerator hous i ng (22). 

After installation of the regenerator (18), a sintered stack of stain- 
less steel wire mesh, each assembly is mounted on the front crankcase. 

The piston (21) is a light weight welded assembly consisting of a base, 
which houses the dry running piston rings and bearing, welded to the piston 
rod and dome. 

The oil scraper/capseal assembly (25) separates each cycle from the 
crankcase. The capseal separates the fluctuating cycle pressure from the 
mean crankcase pressure. The oil scraper is a babbit sleeve radially pressed 
against the rod. There is no pressure difference across the scraper. Its 
only function is to prevent oil from penetrating into the cycle. Both the 
oil scraper and the capseal are mounted in a compliant housing to 
accomodate the slightly eccentric motion of the piston rod due to the 
running clearance of the cross head. 

The cross heads (26) are made from a high strength aluminum alloy. 

The extended length reduces side forces and bridge deflection. The running 
surfaces are treated with a material that is compatible with the Reynolds 
390 bore. 

The life of the engine is determined by the creep failure life of the 
hot parts, and is well in excess of 50,000 hours at full load. 

Table 1 summarizes some of the important features and parameters 



Arrangement: 


Bore; 

Maximum Stroke: 

Overall Length: 

Cross-sec t i ona 1 Dimensions: 
Total Estimated Weight: 
Working Fluid: 

Mean Cycle Pressure: 

Heater Temperature: 

Power Control : 

Heat Transport: 

Gas Containment: 

Oil Containment: 

Materials: 


Four double-acting cylinders symmetrically 
arranged about a common axis. One heat 
exchanger assembly per cylinder. 

56 mm 

48 mm 

635 mm 

Largest cross section is 300 mm in diameter 

75 kg 

Helium 

11 MPa 

800 °C 

Piston stroke variation by means of a 
variable swashplate with a maximum angle 
of 22° 

Sodium heat pipe 

Crankcase pressurized to mean cycle pressure 
and sealed with a rotating shaft seal 

Reciprocating rod oil scrapers 

Iron-base CRM-6D, CG-27 heater tubes 


Table 1 


Important Features and Parameters of the STM4-120 




PERFORMANCE 


A set of engine parameters governing the performance of the 
thermodynamic section was selected for the STM4-120 through an extensive 
and painstaking effort to tailor the engine performance to the design 
approach described above. 

This approach required a high level of efficiency to prevail over 
a wide range of operating conditions, making the engine suitable for 
any duty cycle. Hundreds of different combinations of engine parameters 
were simulated before the optimal set was selected and established as 
the design base of the engine. 

The result is shown, in Figures 12a and 12b, as two performance maps 
at mean cycle pressure levels of 11 MPa and 6.3 MPa respectively. Figure 12a 
shows that between the power levels of 8 kU and 40 kW the shaft efficiency 
is between 45% and 47% (excluding auxiliaries and collector efficiency). 

Over a wide power range the efficiency hardly changes with engine speed. 

For applications requiring less power, a smaller charge of 
helium may be used with very little effect on the efficiency. This is 
shown in Figure 12b for an engine charged to 6.3 MPa to provide power 
output of no more than 25 kW. This will greatly enhance the life of the 
engine. 

Major contributions to the high efficiency of the engine are derived 
from the heat-exchanger-stack configuration, the variable stroke power 
control and the design of the drive mechanism. 

The heat-exchanger-stack configuration adds more freedom to the set 
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List of Symbols 

A Heat transfer area [m*] (for gap: dA = 2wdx) 

g Gap width [m3 

h Heat transfer coefficient [ ] 

m x r 

k Thermal conductivity ] 

m Mass flow [kg/sec] per gap 

Nu Nusselt number (for gap: Nu * hg/k) 

Q Heat flow (Watts) 

T Temperature (°C) 

t Preheater plate wall thickness (m) • 

w Depth of heat exchanger (or preheater) gap (m) 
x Coordinate in the flue gas flow direction (m) 

cp Specific enthalpy (3/kg) 

Subscripts 

f denotes properties of the flue gas 
a denotes properties of the incoming air 
w denotes properties of the heat exchanger wall material 


\ 
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Governing Equations for Preheater Heat Transfer Analysis 

The analysis is performed on a counterflow, fiat plate recuperator with 
multiple parallel gaps of width (g), depth (w) and length ( Jt ). 

Laminar flow is assumed throughout the gaps and hence the Nusselt 
number (based on the gap width) is constant: 

Nus £ g = 3.8 (1) 

The analysis, performed on a single gap, takes into account the varia- 
tion Of fluid properties with temperature, but not with pressure. 

The functional relationship of the thermal conductivity (k) and the 
enthalpy (<p) to the temperature T for both flue gas and incoming air is 
known for the particular fuel used and the equivalence ratio. The latter 
relation takes into account molecular dissociation. 

The energy balance for the flue gas requires that heat loss from the 
flue gas by convection and conduction to the incoming air stream be equal to 
the reduction in enthalpy of the flue gas stream: 

-<T - T ) dA 

<«, ■ i J — nr- - < T f> (2 > 

h x + k... + h 


Observing that dA = 2wdx, substituting (1) in (2) and rearranging, we 
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Similarly, energy balance for the air .stream requires that: 


0Q a . m a d » a <T a ) 


dQ a = -dQj = mjd<pj(Tj) 


hence: 


. <ff a <T a ) . d» f <T f > 
m a dx * ‘"'f dx 


Equations (3) and (6) form a set of two coupled, first order, non- 
linear ODE's for the two unknown functions T^(x) and T a <x). The boundary 
conditions are the known inlet temperatures T.(O) and T (&) of the flue gas 

I A 

and air streams respectively. 

Equations for the Heat Exchaneer/Evaporator Analysis 


The analysis is performed on a heat exchanger with flat parallel gaps 
(of width g, depth w, and length i ) through which flue gas flows laminarly 
at a constant Nusselt number Nu - 3.8. The gas yields its heat to the gap 
walls which are at constant and uniform temperature t . 



As before, energy balance on the flue gas requires: 


dQ { = -h f (T f - T w ) dA = r^ f d<p f (T f ) 


Substitution of 


and of 


in (7) yields: 


Nu 



dA - 2wdx 


m f dx 


k f (T c ) 

* 2w Nu-^-£- (T f - T w ) 


= 0 


This first order, non-linear ODE in the unknown function T^(x) forms the 
mathematical model of the evaporator/heat exchanger. The boundary condition 
is the known inlet temperature T f (o) which is the flame temperature. 
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EHS Analysis Code Listing 


Fuel Data 


atomic ratio of hydrogen to carbon in fuel 
specific lower heat value of fuel, cal./gr, 


Operating Data 


MFUEL 


equivalence ratio 

atomizing air to intake air ratio 

initial guess for preheated air temperature 

fuel consumption rate, gr./sec. 

outside wall temperature of evaporator °K 

air intake temperature °K 


Geometric Data; Heat Exchanger (see figure) 

EXL length, m 

WEX width, m 

NGEX number of flue gas passages ("gaps") 

GEX gap width, m 

NUEX Nussel number 


wall thickness, m 

thermal conductivity of the metal, w/m °K 


Geometric Data: Preheater (see figure 


PL length, m 

WP width, m 

NGP number of flow passages for each stream 

GPA gap width, air side, m 

GPF gap width, flue side, m 

NUP Nussel t number 


Cold Start Data 


density of the metal, gr/cm^ 


specific heat of the metal, cal./gr. °K 


effective gauze thickness: the thickness of a solid sheet 
of metal of the same weight per unit area as the gauze 
stack, m 






Numerical Data 


NEX number of discrete steps in heat exchanger calculation 

NPA number of discrete steps in preheater calculation 

EPSP maximum convergence error in preheater iterative 

procedure, relative 

EPS maximum convergence error in EHS iterative procedure 

relative 
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SAMPLE RUNS - AUGUST 8, 1985 


dun. # l (o} <k * c y£’ fed** 4*yte 646*€***t *n G$fO* &** 


1 . 8040 

10179. 


1 . 2500 

. 40000E-01 

949. 90 

. 1 2000 

. 60000E-01 

61 . 

. 18000 

. 60000E-0 1 

78 . 

8.0000 

. 12000 

-. moo. 

100 100 

. 20000E-02 . 

20000E 


2.0000 1075.0 301.00 

>E-02 3.8000 . 30000E-03 21.600 

>E-02 . 13000E-02 3.8000 


PREHEATED AIR TEMP. 

949.9 

K 

FLAME TEMPERATURE 

2374. 9 

K 

HEAT EXCHANGER OUTPUT TEMP. 

1096. 3 

K 

HEAT OUTPUT 

74.661 

kw 

EXHAUST OUTLET TEMP. 

569.9 

K 

E.H.S. EFFICIENCY 

87.675 

7. 

COLD START PENALTY 

45. 389 

gr 

PRESSURE DROP 

2. 193 

KPa 

FLOW FRICTION POWER 

296.791 

w 


623.38 


841.16 


1388. 4 


B>A SPRINT EHS.FOR 


Not ready error reading drive A 
Abort, Retry, Ignore? R 

BsEHS .FOR is currently being printed 



V 



B> 


3t 0.12,0.06, 61 1 . 3E-3, 3.8,0. 3E-3, 21.6, 

4s 0.18,0.06, 78 1 . 3e-3, 1 . 3e-3, 3. 8, 

5s* 8.0,0. 12,0. 537E-3, 

6s 100 100 0.002,0.002 

*E 

B>EHS 

1.8040 10179. 

1 . 25«.".» . 40000E-01 949.90 2.0000 1075.0 301.00 


.12000 . 60000E-01 61 . 13000E-02 3.8000 

.18000 . 60000E-01 78 . 13000E-02 . 13000E-02 


8 . 0000 . 1 2000 . 53700E-03 

100 100 . 20000E-02 . 20000E-0i? ' 


PREHEATED AIR TEMP. 

FLAME TEMPERATURE 

HEAT EXCHANGER OUTPUT TEMP. 

HEAT OUTPUT 

EXHAUST OUTLET TEMP. 

E.H.S. EFFICIENCY 
COLD START PENALTY 
PRESSURE DROP 
FLOW FRICTION POWER 
623.38 841.16 1388.4 


949.9 K 
2374.9 K 
1096.3 K 
74.661 kw 

569.9 K 
87.675 7. 
87.361 gr 

2.193 KPa 
296.791 w 


. 30000E-03 2 1 . 600 
3 . 8000 
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1 VJ->- 


End of input file 

*6 

6:* 100 100 0.002,0.002 
6s* 200 200 0.002,0.00^ 


B >EHS At 
1.8040 
1 . 2500 
. 12000 
. 18000 
8.0000 
200 200 


Ret*. & 3 Atcf 

10179. V 


. 40000E-01 949.90 2.0000 1075.0 301.00 

. 60000E-01 61 . 13000E-02 3.8000 . 30000E-03 21.600 

. 60000E-01 78 . 13000E-02 . 13000E-02 3.8000 

. 1 2000 . 53700E-03 


.18000 . 60000E-01 78 .130 

8.0000 .12000 . 53700E-0 

200 20C{ . 20000E-02 . 20000E-02 

PREHEATED AIR TEMP. 

FLAME TEMPERATURE 

HEAT EXCHANGER OUTPUT TEMP. 

HEAT OUTPUT 

EXHAUST OUTLET TEMP. 

E.H.S. EFFICIENCY 
COLD START PENALTY 
PRESSURE DROP 
FLOW FRICTION POWER 

624 . 00 837 .51 1 383 . 4 


950.5 K 
2374.9 K 
1097.2 K 
74.613 kw 
570.4 K 
87.618 7. 
87. 187 gr 
2.179 KPa 
293.849 w 


3:* 0.12,0.06, 61 1 . E-3, 3. 8, O. 3E-3, 21 . 6, 

4:* 0.18,0.06, 78 1 . 3e-3, 1 . 3e-3, 3. 8, 

4s* 0.18,0.06, 78 1 . e-3, 1 . e-3, 3. 8, 


B >EHS 
1 . 8040 


ztH fty 


10179. 


1 . 2500 . 40000E-0 1 949 . 90 

.12000 . 60000E-01 61 . 1 00' 

.18000 . 60000E— 01 78 . 100 

8.0000 .12000 . 53700E-0 

200 200 . 20000E-02 . 20000E-02 

PREHEATED AIR TEMP. 

FLAME TEMPERATURE 

HEAT EXCHANGER OUTPUT TEMP. 

HEAT OUTPUT 

EXHAUST OUTLET TEMP. 

E.H.S. EFFICIENCY 
COLD START PENALTY 
PRESSURE DROP 
FLOW FRICTION POWER- 
645. 59 818.03 1319.2 


49.90 2.0000 1075.0 301.00 

61 . 1 0QQ0E-02 3.8000 . 30000E-03 21.600 

78 . 1 0000E-D2 . 10000E— 02 3.8000 
53700E-03 
0000E-02 

976.9 K 
2382.0 K. 

TEMP. 1082.8 K. 

76.464 kw 
532.0 K 
89.792 7. 

85.233 gr 
4. 493 KPa 
578.828 w 


Xa£a A /. «•- 
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m 

b 

Q 

A 

v 

P 

X 

P 

g 

T 

r 

Re 

f 

0 

e 

K 

F 

D 

*,y 
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Subscripts 

Jl 

v 

w 


List of Symbols and Abbreviations 

Mass flow rate 

Width 

Length 

Heat flow 

Area 

Velocity 

Density 

Heat of vaporization 
Pressure 

Acceleration of gravity 

Shear stress 

Radius 

Reynolds number 
Drag coefficient 

Angle of liquid flow to horizontal 

Porosity 

Permeability 

Friction coefficient 

Dynamic drag coefficient 

Position 

Profile coefficient for momentum flow 


Liquid 

Vapor 

Wick 

Hydraulic 
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Define friction coefficient, F^; and substituting into Equation 13; 


F * M * 
1 ' 


(U) 



-F^Q ( 1 - -y- ) - Pgg sin 0 


(15) 


Integrating Equation 15 yields an expression for the liquid pressure drop 
(including the hydrostatic term) as required for evaluating the capillary pressure in 
Equation 18 in the text. 

From conservation of momentum, the vapor pressure drop can be expressed 
as follows: 


dP 
v 

dy 


-(f Re ) p m 
v v r v v 


2A r 7 p 
v v v 


- - 0 


2m 

v 

%# n % 


dm 
v 

dy 


(16) 


The construction of the fin is such that the vapor path is packed with a porous 
screen structure to provide structural support to the fin walls. This support structure 
is analagous to the screen wick used for the liquid flow, with rega; d to the flow 
resistance since it can be shown that the vapor flow is laminar and incompressible 
(see Appendix D). Therefore, by introducing the vapor screen permeability, K y , and 



vapor screen porosity e v , Equation 16 becomes; 


dP 


= - [ 


] "k* Q * [ a rz i i ] 2 (“u7") Q 07) 


dy ' ' L KA. X p.. J b v L V j 


v v 


B 

n 

V V 


Introducing the vapor friction coefficient, F y and the vapor drag coefficient, D y yields 


= -F y Q ( f) . 2D„ -§- 2 <V> 


(18) 


v b ' b 

Integrating Equation 18 yields an expression for the vapor pressure drop as required 
for evaluating the capillary pressure in Equation 18 in the text. 
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PRELIMINARY DESIGN DRAWINGS 


rlt 






PM Hf AT If *'<bTH 

ift.i —(t*n 








gggl 



y>-r 


















